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ABSTRACT 

Aims. We present new high resolution (R > 50,000) absorption measurements of the Nal doublet (5889 - 5895A) along 482 nearby sight-lines, 
in addition to 807 new measurements of the Call K (3933A) absorption line. We have combined these new data with previously reported 
measurements to produce a catalog of absorptions towards a total of 1857 early-type stars located within 800pc of the Sun. Using these data we 

, have determined the approximate 3-dimensional spatial distribution of neutral and partly ionized interstellar gas density within a distance-cube 

■ of 300pc from the Sun. 

' Methods. All newly recorded spectra were analyzed by means of a multi-component line profile-fitting program, in most cases using 
simultaneous fits to the line doublets. Normalized absorption profiles were fitted by varying the velocity, doppler width and column density for 

' all intervening interstellar clouds. The resulting total column densities were then used in conjunction with the Hipparcos distances of the target 
stars to construct inversion maps of the 3-D spatial density distribution of the Nal and Call bearing gas. 

Results. A plot of the equivalent width of Nal versus distance reveals a wall of neutral gas at ~80pc that can be associated with the boundary 
wall to the central rarefied Local Cavity region. In contrast, a similar plot for the equivalent width of Call shows no sharply increasing 
absorption at 80pc, but instead we observe a slowly increasing value of Call equivalent width with increasing sight-line distance sampled. 



I Key words. Galaxy:solar neighborhood - ISM:atoms - ISM: clouds 



1. Introduction 



Knowledge of the spatial distribution, dynamics and the asso- 
ciated physical and chemical state of interstellar gas can pro- 
vide important insights into the subtle interplay between the 
evolution of stars and their exchange of material to and from 
the ambient interstellar medium (ISM). Although much valu- 
able work in this field has been achieved through numerous 
21cm radio and infrared surveys of our Galaxy, knowledge of 
the physical state and spatial distribution of of the interstellar 
gas within 300pc of the Sun is still far from complete. The low 
density region surrounding our Sun to ^ lOOpc in all direc- 
tions is termed 'the Local Cavity' (LC) and has been shown to 
be largely free of cold and dense gas, although it does contain 
many partially ionized diffuse cloudlets often referred to as 'lo- 
cal fluff (Redfield & Linskv l2008] l. The actual physical state 
of the gas that presumably fills the remaining volume in the 
Local Cavity is still much debated. The prevailing view is that 
it is mainly filled with a very low density, highly ionized (mil- 
lion K) gas, as inferred by observations of the diffuse soft X- 
ray (0.25keV) background emission (Snowden et al. 1998) and 
also by the detection of local absorption from the high ioniza- 



tion OVI (1032A) line by Savage & Lehner ( [2006l l. However, 
several authors have raised some serious arguments against the 
prevailing interpretation of these data, such that the LC gas 
may be of a far lower temperature or in a state far from ioniza- 
tion equilibrium (Breitschwerdt |200Tl Koutroumpa et al. 120071 
Barstow et al. 2008 ). Based on the demonstration that most of 
the unabsorbed soft X-ray brightness at low galactic latitudes 
is actually solar wind charge-transfer X-ray emission gener- 
ated within the heliosphere, Koutroumpa et al. ( 120091) have sug- 
gested that the LC may contain far cooler or more tenuous gas 
at low altitudes (| z \ < 60 - lOOpc) with the soft X-ray emis- 
sion seen from | & | > 30° arising from an 'external' hot gas 
and from the overlying halo (Welsh & Shelton |2009l l. 

Our unique placement within this low density interstellar 
cavity allows us to directly measure its spatial extent and phys- 
ical state through high spectral resolution absorption measure- 
ments that are mostly uncontaminated by intervening dense 
line-of-sight interstellar features. Over the past 20h- years we 
have been carrying out high spectral resolution (R ~ 100,000) 
absorption observations of the local neutral and partially ion- 
ized interstellar gas using both the Nal D-line doublet at 5890A 
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and the Call-K line at 3933A (Lallement et al. [19861 Vallerga 
et al. [T993l Welsh et al. [19941 Sfeir et al. [T999l Lallement et al. 
I2003I I. These, and many other, observations have allowed us to 
construct a picture of the distribution of neutral gas in the local 
ISM in which the LC is connected though interstellar tunnels to 
other surrounding cavities, as predicted in the model of the ISM 
by Cox & Smith (1974). In particular, maps of the distribution 
of Nal absorption have revealed 3 large interstellar features that 
dominate the appearance of the local ISM within 200pc. These 
are: (i) the 50pc diameter and 200pc long extension to the LC in 
the direction of the star /3 CMa (Welsh|T99TJ, (ii) the extension 
of the rarefied LC into the lower galactic halo to form an open- 
ended Local Chimney feature (Welsh et al. 119991) . and (iii) the 
apparent connection of the LC with the large Loop I superbub- 
ble at a distance of 80pc (Welsh & Lallement 2005). We 
attribute the name "Loop 1 superbubble" to the cavity centered 
at longitude 345 deg, although the link between this cavity and 
the radio feature is far from well understood. It seems likely 
that the existence of these 3 features is somehow linked to the 
formation history of the LC itself. 

The need for accurate maps of the 3-D density distribution 
of neutral local gas has been demonstrated in wide variety of 
recent studies (Oegerle et al. 2005; Smith et al. 120051 Combet 
et al. [20051 Posselt et al. i2007.: Wolleben (20071). Such maps 
are particularly relevant for studying the interaction regions 
between local gas residing in the galactic plane, gas in adja- 
cent interstellar cavities and local gas being ejected into, and 
falling from, the lower galactic halo. It is important to note that 
only when absorption data are combined with parallax mea- 
surements can the distances to gas clouds be accurately derived. 
Maps of the distribution of interstellar gas derived from HI, 
CO and the IR are of two-dimensions and cannot provide ac- 
curate distance estimates to gas clouds. Even so, the accuracy 
and spatial extent of our existing absorption maps are presently 
limited by (a) the number of sight-line directions sampled, (b) 
the sampling length along and between each sight-line and (c) 
the distance over which each sight-line extends. At present the 
most accurate maps of the density distribution of local neutral 
gas are those presented by Lallement et al. (120031). Some 1005 
sight-lines within ~ 250pc were observed using Nal absorp- 
tion in order to produce 3-D neutral gas density distribution 
maps to a distance of ^ 150pc in most directions. The galac- 
tic distribution of these 1005 sight-lines is shown in Figure 1 of 
Lallement et al. (I2003I I. hereafter Paper 1 . Many of these obser- 
vations traced the extent of Gould's Belt, since only early-type 
stars were used as background absorption sources. We note that 
in Paper 1 the typical global sampling area on the sky was ^ 1 
target per 50 sq. deg., with significantly less coverage at high 
galactic latitudes. 

No equivalent maps of the density distribution of partially 
ionized gas, as traced by the interstellar Call absorption line, 
currently exist for the local ISM. This is particularly important 
for sight-lines ^ SOpc, since the LC gas is thought to mainly ex- 
ist in an ionization state higher than that which can be probed 
by the Nal ion (Welsh et al. 1 19941 ). Many diffuse and partly 
ionized cloudlets are known to exist within this volume and 
previous studies have demonstrated that the Call K-line is an 
excellent tracer of absorption of this physical phase of the local 



gas (Lallement, Vidal-Madjar & Ferlet ll986l Crawford, Craig 
& Welsh [T997I) . In fact, in a comparison study of high resolu- 
tion ground-based Call K-line and UV (Fell and Mgll) absorp- 
tion line data towards stars within lOOpc by Redfield & Linsky 
(120021) . it was found that most (but not all) sight-line compo- 
nents were observed in both the Call and UV lines to within 
a velocity error of only ±1.5 km s^^. This evidence suggests 
that the Call line is a good indicator of the absorption structure 
that one might expect from the warm and ionized gas sampled 
by the UV lines. The good agreement between these absorp- 
tion velocities clearly suggests a co-location for the formation 
of both the optical and UV lines within these same warm and 
diffuse local gas clouds. Furthermore, UV observations of stars 
within the (3 CMa interstellar tunnel suggest that although the 
neutral gas density along this 200pc long interstellar feature 
may be very low, this is not the case for ionized gas components 
(Dupin & Grv 11998]) . Clearly it would be instructive to com- 
pare maps of the spatial distribution of Call absorption (which 
is sensitive to partially ionized gas) and those derived from Nal 
observations (which trace cold and neutral gas) for the local 
ISM. 

In this Paper we report on the results of an extended survey 
of Nal and Call absorption lines recorded at high spectral res- 
olution (R > 50,000) towards many newly observed early-type 
stars mostly located within 800pc of the Sun. We present new 
absorption measurements of the Nal doublet recorded along 
482 sight-lines, in addition to 807 new measurements of Call 
absorption along 807 sight-lines. When added to similar Nal 
and Call absorption measurements reported in the literature, 
we are now able to present new maps to 300pc of the 3-D 
distribution of neutral and partially ionized gas towards 1678 
stars (at Nal) and 1267 stars (at Call). Although these new data 
confirm many of the findings of Paper 1, they also show sev- 
eral new interstellar gas features now revealed through the in- 
creased level of sampling of local ISM sight-lines, especially at 
high galactic latitudes and in regions with distances in the 150 
- 350pc regime. The 3-D maps of local Call absorption are the 
first of their kind to be presented, and in several respects they 
mimic many of the interstellar features seen in the equivalent 
maps of local Nal absorption. However, there are a significant 
number of regions where Call absorption is high and Nal ab- 
sorption is low, and vice-versa. Such variations in the observed 
column density ratio of N(NaI)/N(CaII) are influenced by the 
conditions in the local stellar ionization field. Values of this ra- 
tio for stars near to the galactic plane in each of the 4 galactic 
quadrants are also presented. Finally, as a by-product to de- 
termining the (Nal and Call) gas column densities along each 
sight-line we have also determined the absorption velocity (i.e. 
cloud component) structure towards each target. These results 
will be published in a separate Paper that deals with the kine- 
matics of local gas clouds within ^ 150pc (Lallement et al. 
[20T0I) . 

2. Observations and Data Reduction 

We present observations of both the interstellar Nal D-line 
doublet at ~ 5890A and the Call K-hne at 3393A recorded 
in the absorption spectra of early-type stars with distances < 
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0.8 kpc. These data were obtained during observing runs per- 
formed over the period 2003 - 2008 using the following in- 
strumentation (i) the Aurelie spectrograph at the 1 .52 m tele- 
scope of the Observatoire de Haute Provence (France), (ii) 
the Hamilton echelle spectrograph on the 0.9m coude feed 
telescope at the Lick Observatory (USA), (iii) the Hercules 
echelle spectrograph on the 1.0m McLellan telescope at the 
Mt. John Observatory (New Zealand), (iv) the FEROS coude 
echelle spectrograph on the 2.2m telescope at the European 
Southern Observatory (Chile) and (v) the GIRAFFE fiber-fed 
echelle spectrograph at the 1.9m Radcliffe telescope of the 
South African Astronomical Observatory (RSA). The spec- 
tral resolution for these interstellar data was 3 km s^^ for the 
Observatoire de Haute Provence (OHP) data, 5 km s^^ for the 
Lick Observatory data, 4 km s^^ for the Mt. John Observatory 
(MJO) data, 6 km s^^ for the FEROS ESO data and 7.5 km s^^ 
for the South African Astronomical Observatory (SAAO) data. 

Interstellar sight-lines were generally selected using the cri- 
teria that the background stellar targets possess: (i) a Hipparcos 
catalog distance (ESA U 9971 ) normally less than 400pc and with 
an associated relative standard parallax error smaller than 0.3 
(however several targets with distances up to 800pc were also 
included for certain sight-lines of specific interest), (ii) a spec- 
tral type earlier than A5 V, (iii) a stellar rotational velocity > 25 
km s~^, (iv) a galactic position that filled in areas of the sky 
not well sampled by the targets listed in Paper 1, and (v) a suf- 
ficiently bright visual magnitude that a well-exposed spectrum 
could be gained within a 60 minute exposure. We also avoided 
observations of known spectroscopic binary stars, since inter- 
stellar features are often difficult to identify within their stellar 
absorption profiles. In Table 1, available on the CDS, we list in- 
formation on the newly observed targets that includes their HD 
number, galactic longitude and latitude, Hipparcos distance (in 
pc) and the observatory (OHP, Lick, MJO, ESO and SAAO) 
where the data was taken, together with the year of observa- 
tion. 

All of the data, except for the MJO and ESO observations, 
were reduced using software routines outlined in detail in Sfeir 
et al. (I1999I I. These routines involve: (i) division of the raw 
data by a flat-field, (ii) removal of instrumental background 
light and cosmic rays, (iii) removal of telluric water vapor lines 
(which particularly affect the Nal D2-line at A 5890A ) us- 
ing either a synthetic telluric transmission spectrum as detailed 
in Lallement et al. (119931 1 or by division by a purely stellar- 
atmospheric absorption spectrum of an unreddened (nearby) 
B star observed at a similar altitude and time, and (iv) spec- 
tral order extraction and wavelength calibration using Th-Ar 
lamp spectra recorded several times throughout each observing 
night. Similar data reduction processes were performed on the 
Mt. John Observatory and ESO data using their own in-house 
standardized echelle data extraction software packages. 

The wavelength calibrated spectral data were subsequently 
fit with a high order polynomial in order to define a local stellar 
continuum placement, such that the equivalent widths (Wa) of 
the Nal and Call absorption lines could be measured from these 
resultant intensity profiles. The errors on the measured values 
of W\ are dominated by the statistical uncertainty on the mea- 
sured counts within each line profile and the systematic uncer- 



tainty in the local continuum placement (see Welsh et al. 119901 
for a fuller discussion). These two sources of error are added 
in quadrature to derive the measurement error Typically, for a 
well-exposed spectrum, a measurement error of ~ 5% is nor- 
mal for our data. Upper limits to the value of equivalent width 
where no absorption line was detected with significance were 
derived from conservative estimates of the strength of a poten- 
tial absorption feature appearing at a level > 2.5-(j above the 
rms value of the local continuum. For the majority of our well- 
exposed spectra (i.e. S/N > 30:1) this typically resulted upper 
limit values of Wa < 2.5mA for both the Nal-D2 and Call-K 
lines. The measured values of Wa for the Nal (D2 & Dl) and 
the Call-K lines are listed in Table 1 for our presently observed 
sight-lines. 

For the majority of targets with distances 300pc both the 
interstellar Nal and Call absorption lines are not greatly satu- 
rated. Thus, we were able to accurately fit these profiles with 
absorption components (i.e. gas clouds) using the line-fitting 
procedure described in Sfeir et al. (I1999I I. A best-fit theoret- 
ical absorption profile is generated that is characterized by 3 
parameters: (i) an interstellar cloud (heliocentric) velocity, V 
(ii) a doppler velocity dispersion parameter (6-value) and (iii) 
a cloud component column density, N(Nal) or N(Call). The 
best-fit models were further constrained by fitting both of the 
Nal D-lines simultaneously (with the fit being biased towards 
the stronger D2 line). For the ESO data, in which both the Call 
K and H-lines were available, the best-fit for the doublet was 
similarly constrained. 

As an indication of the quality of the data recorded with 
the spectrograph-Htelescope systems at each of the 5 observato- 
ries, in Figures 1- 5 we show some typical model fits to the ob- 
served interstellar Nal and Call absorption profiles recorded to- 
wards stars with distances in the 50pc - 500pc range. The best- 
fit model parameters (i.e. V, b and N(Nal)) and/or N(Call) for 
the (several) cloud components required to fit each of these ab- 
sorption profiles are listed in a box beneath each of the profiles 
in these Figures. The fits were performed using the minimum 
number of absorption components, with the addition of extra 
components only being deemed necessary until the residual 
between the model fit and the data points was approximately 
equal to the rms error of the data. In general (dependent on the 
spectrograph used), the accuracy of the placement of each of 
the cloud component velocities within these line-fits is ~ ±0.5 
km s~^. 

The highest resolution data (see Figure 1) was obtained for 
the Call K-line at the OHP, and hence it is no surprise that 
these data generally require more absorption components to 
fit the observed profiles. In contrast, the data gained with the 
lower resolution FEROS system on the 2.2m telescope (shown 
in Figure 4) is of the highest S/N ratio and allows access to both 
the Nal and Call doublets, such that the model fits can gener- 
ally be constrained to a better degree of confidence. However, 
for the purposes of our present interstellar study we only re- 
quire use of the total equivalent width, Wa(D2) and Wa(K), 
and the summed total column density values, N(NaI)tot and 
N(CaII)tot, measured towards each star. These values are listed 
in Table 1 . The best-fit values of Nal and Call column density 
for each cloud component for aU of the newly observed targets 
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will be presented in a forthcoming paper which will concen- 
trate of the velocity structure of the local ISM within 150pc 
(Lallement et al |20T0] l. 

In order to increase the spatial sampling of interstellar gas 
absorption in the local ISM, in addition to the newly obtained 
data listed in Table 1 we have also included Nal and Call inter- 
stellar measurements (of equivalent width and associated total 
column density) towards sight-lines < 800pc that have been 
previously published in the literature. The corresponding data 
for each of these observations are also listed in Table 1 along 
with the appropriate reference for the particular observation. 
For cases in which only values of equivalent width are quoted 
in the literature, we have derived a corresponding value of total 
column density (reported in Table 1) using empirical curves of 
growth for both Nal and Call derived from all of the other data. 
We have not included interstellar data for targets with known 
(or suspected) circumstellar gas, since this can provide a sig- 
nificant (and variable) contribution to the values of Nal and 
Call column density (Welsh et al. I1998I I. One important case 
of this type is the star HD 225132 (2 Cet), which appeared in 
the Nal maps of Paper 1 as a nearby (d = 70pc) condensation of 
cold gas lying within the confines of the LC. Recent high res- 
olution measurements of Nal and Call absorption by Redfield 
et al. (I2007I I. which have subsequently been confirmed by our 
own measurements, revealed a far smaller value of Nal column 
density for this sight-line compared to that reported in Welsh 
et al. (|1994| l. In addition, our own Call measurements reveal a 
component at V = +20 km s^^ not seen by Redfield et al. We 
have therefore removed this star from our compilation in Table 
1. In total we now present 1678 Nal and 1267 Call interstellar 
sight-line absorption measurements, 1088 of which have com- 
mon sight-line measurements for both of these ions. 

In Figures 6(a) and 6(b) we show the galactic distribu- 
tion of the respective 1678 Nal and 1267 Call sight-lines sam- 
pled. Although the sky-coverage is significantly increased since 
Paper 1, we note some areas of the sky (e.g. I = 295°, b = -45° 
and I = 160°, b = -75°) remain sparsely sampled. This is pri- 
marily due to the paucity of suitably bright O, B and A-type 
stellar targets in these directions that can be observed in a rea- 
sonable exposure time with the presently available instrumen- 
tation. Therefore any conclusions derived from the mappings 
of interstellar absorption towards both of these galactic direc- 
tions must be viewed with some caution. We also note that our 
selection of suitable targets is biased against heavily reddened 
and visually faint (m^, ^ 6.5) sight-lines, due to the practical 
restriction of ~ one hour of exposure time per recorded spec- 
trum. 

3. Interstellar Absorption Mapping 

The measurements of Nal and Call equivalent width and col- 
umn density measured along the 1857 sight-lines allows us to 
probe the absorption characteristics of both neutral and par- 
tially ionized local interstellar gas. Our underlying assumption 
is that the Nal ion is a good tracer of of the total amount of neu- 
tral and cold (T < lOOOK) interstellar gas along most galactic 
sight-lines (Hobbs ,1978) . Similarly, the Call ion is assumed 
to trace both neutral and warmer (T < 10,000K) partially ion- 



ized gas clouds. Previous high resolution studies of Nal and 
Call absorption suggest that the latter ion has larger absorp- 
tion line-widths for components seen in Nal absorption at the 
same velocity (Welty, Morton & Hobbs [19961 1. This has been 
interpreted as the Call gas components occupying a larger vol- 
ume than those of Nal, with an associated higher temperature 
and/or more turbulent velocity than the volume of gas occupied 
by Nal. In our present study, in which we average the absorp- 
tion properties of gas along many interstellar sight-lines, we 
assume that the derived spatial distributions of both Nal and 
Call are directly comparable. We note that Call absorption not 
only arises in cold (and relatively dense) gas in which Ca may 
be heavily depleted onto dust grains and is the dominant ion- 
ization state, but it can also be present in warmer and lower 
density gas clouds where Ca is less depleted and in which Call 
is a trace ionization state. 

Our primary aim in this paper is to determine the 3-D spa- 
tial distribution of interstellar gas (as characterized by Nal 
and Call absorption) to a distance of ^ 300pc from the Sun. 
Following our previous work on this subject (Sfeir et al l 19991 
Lallement et al. '2003), we produce plots of the spatial distribu- 
tion of gas based on the inversion of column densities using a 
method originally devised by Vergely et al j2001l l. This method 
of plotting, together with some of its advantages and short- 
comings, has been described fully in Lallement et al. (I2003I I. 
However, before presenting such gas density maps we firstly 
discuss the measurements of the equivalent widths of the inter- 
stellar Nal and Call lines. 

3. 1. Nal and Call Equivalent Width Values 

In Figure 7 we show a plot of distance versus the total equiv- 
alent width of the Nal D2-line, Wa(D2), for targets with dis- 
tances < 400pc. Note that this plot does not include sight-lines 
that have reported values of N(Nal), but no associated informa- 
tion on the measured equivalent width. This figure, in agree- 
ment with a similar plot shown in Sfeir et al. (119991 ). shows 
that very little measurable Nal absorption can be detected for 
distances (in all galactic directions) up to ^ 80pc from the Sun. 
This general lack of appreciable Nal absorption (Wa(D2) 5 
mA) provides clear evidence that the volume of the LC is es- 
sentially free of major condensations of cold and dense neu- 
tral interstellar gas. Beyond 80pc the level of Nal absorp- 
tion rises sharply over a short distance of 20pc to reveal the 
presence of a dense 'wall' of neutral gas (Wa(D2) > 20 mA) 
that surrounds the LC in many galactic directions (Welsh et 
al. 1 19941 1. Beyond the LC neutral boundary, which is normally 
characterized by a single narrow velocity component in the Nal 
absorption profiles, more and more Nal absorbing (cold) neu- 
tral interstellar clouds are encountered as sight-line distances 
increase. In several galactic directions with distances > 300pc 
the Nal D-lines become increasingly saturated and we therefore 
generally see a slower increase in equivalent width value with 
greater distance sampled. However, we also note the presence 
of several sight-lines with distances > 250pc that have minimal 
associated Nal absorption (plotted as open circles in Figure 7). 
These are sight-lines which extend into the lower galactic halo 
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through the openings of the Local Chimney, a region of known 
low neutral gas density (Crawford et al |20021 l. The nearest stars 
with distances less than 65pc, but with 'anomalously high' val- 
ues of Wa(D2) > 10 mA , are HD 184006 (d = 38pc), HD 
159170 (d = 48pc), HD 186882 (d = 52pc), HD 105850 (d = 
56pc), HD 96819 (d = 58pc) and HD 129685 (d =63pc). 

In contrast with studies of interstellar Nal, previous inves- 
tigations of Call absorption have shown that the local ISM is 
composed of many diffuse and partially ionized warm (T 
7000K) cloudlets that possess a complex velocity structure, 
even over distances as short as 5pc (Crawford, Lallement & 
Welsh 1 19981 Lallement et al. 119861) . These cloudlets are ran- 
domly scattered within 50pc of the Sun and although are gen- 
erally not revealed by Nal absorption measurements, they have 
been extensively studied at ultraviolet wavelengths (Redfield & 
Linsky|2008 ). In Figure 8 we show that the behavior of Call ab- 
sorption as a function of distance (to 400pc) does not follow the 
same pattern that was found in Figure 7 for Nal measurements. 
For the majority of sight-lines within ^ lOOpc we mostly find 
WA(Call-K) < 15 mA, with no sharply increasing value of ab- 
sorption at ~ 80pc that might feasibly be associated with the 
presence of gas associated with the boundary wall to the LC. 
Instead there seems to be a wide-spread local distribution of 
partially ionized Call regions with a range of low integrated 
Call absorption strengths that extends to at least lOOpc. Beyond 
^ lOOpc we observe a slowly increasing value of WA(Call-K) 
with the sight-line distance sampled. Even by 200pc there are 
still several sight-lines in which very little (< 3mA) interstellar 
Call is present. At a distance of ~ 300pc a value of WA(Call-K) 

< 100mA is typical for most of the galactic directions sampled, 
although values in excess of 200mA are found for a few higher 
density regions. 

Welsh et al. ( 119971 have presented plots of the total col- 
umn density of Call versus the total column density of neu- 
tral HI which show that both quantities are broadly correlated. 
However, there is significant scatter in this relationship due to 
the fact that along any given interstellar line-of-sight the mea- 
sured amount of Call-bearing gas is closely linked to the ambi- 
ent ionization conditions, the amount of dust depletion and the 
size and number of gas clouds encountered. The Local Cavity 
is an anomalous region of interstellar space that contains many 
randomly scattered warm and partially ionized low column 
density Call-bearing cloudlets (each of W>,(CaII-K) ^ 5 mA, 
Redfield & Linsky |2002l l, probably mostly of similar physical 
size. Based on our fitting of Call absorption profiles for stars 

< lOOpc, the number of 'cloudlets' encountered along a typical 
interstellar sight-line is generally < 4. Thus, the presently mea- 
sured total equivalent width of Call along sight-lines < lOOpc 
is small (< 15 mA) and variable, dependent on the number of 
warm clouds encountered. For cold interstellar clouds that exist 
beyond the rarefied Local Cavity, it is generally assumed that 
the warm Call gas is present in the 'outer-skins' of these clouds 
(Vallerga et al l 1993] Welty et al [T996l ). Thus, although the neu- 
tral wall boundary to the LC may contain significant amounts 
of cold Nal gas compared with sight-lines lying within the LC, 
the corresponding increase in the equivalent width of Call due 
to the outer skin of this boundary wall is far less pronounced. 
Hence, we do not observe a sharply increasing value in Call 



absorption at 80pc due to the presence of the boundary wall 
of the LC. For greater sight-line distances in the general ISM 
(i.e. d>> lOOpc), the measured Call equivalent width should 
be more closely linked to the number of clouds (or cloud en- 
velopes) sampled along a particular sight-line as opposed to 
the actual distance sampled. Thus we expect to observe a grad- 
ual (but variable) increase in Call absorption equivalent width 
with distance, as is demonstrated in the increase of N(Call) 
with N(HI) by Welsh et al (fT997] l. 

We note that there are only 3 sight-lines within 75pc with 
anomalously large values of WA(Call-K) in the 20 - 30mA 
range. These are the sight-lines towards the stars a Oph (d = 
14pc),HD 145570 (d = 5 Ipc), and HD 1 13852 (d =70pc). Also, 
the sight-hnes towards HD 5487 (d = 259pc), HD 85504 (d = 
322pc) and HD 75855 (d = 446pc) have anomalously low lev- 
els of Call K-Une absorption (i.e. log N(Call) < 10.5 cm~^). 
These are all high latitude sight-lines that extend through the 
openings of the low density Local Chimney into the overlying 
galactic halo. 

3.2. The Volume Density Distribution of Nal and 
Call 

Before progressing to a discussion concerning the creation of 
3-D density maps of the spatial distribution of both interstellar 
Nal and Call gas density, it is informative to firstly investigate 
the global variation of volume density of both ions as a function 
of distance sampled. In Figure 9 we show the plot of Nal vol- 
ume density, njva/, (i.e. N(NaI)/distance) versus distance (pc), 
with upper limit values being plotted as open circles and mea- 
sured values plotted as filled circles. We note that the lowest 
values of njva/ are mostly found for distances < 200pc (i.e. 
the region of the ISM dominated by the effect of the low den- 
sity LC), although there are several sight lines (mostly at high 
galactic latitudes) with distances up to 600pc that have simi- 
larly small volume density values. In the 1 10 - 200pc distance 
range there are about a dozen stars with slightly higher than 
average values of npfai in the range 10^^ - 10^^ cm^"^. The 
majority of these stars lie in the longitude range / = 350° - 50°, 
in the general direction of the galactic center For sight-lines 
> 300pc we see a general trend in the data such that values 
of n^ai in the 10^^ to 10^^" cm"'^ range are representative 
for most distances up to 800pc. We interpret these latter data 
points as indicating that the general ISM has a wide range of 
Nal volume density values, each of which is critically depen- 
dent on the number and density of neutral interstellar clouds 
and interstellar cavities encountered over any given sight-line 
distance. 

In Figure 10 we plot the volume density of Call, ncaii, as 
a function of distance, whose values for sight-lines < lOOpc 
show a complex pattern of behavior For distances < 30pc we 
measure both very high and very low values of ncaii, depen- 
dent on whether a particular sight-line encounters any of the 
local fluff clouds. For distances in the 30 - lOOpc range, al- 
though there are still sight-lines with very low values of ncaii 
< 10^^"cm^^, a value of ~ 5 x lO'^^cm^^ seems fairly rep- 
resentative for most interstellar directions. However, for dis- 
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tances > lOOpc higher values of ncaii are generally encoun- 
tered in the ISM, such that a volume density of '--^10~^cm~^ 
is typical for most sight-lines with distances up to 800pc. This 
would suggest that beyond ^ lOOpc the density distribution of 
Call bearing gas is fairly uniform throughout the ISM. This is 
in accord with the findings of Welsh et al. ( 1 19971 1 and Megier 
et al. (I2005I I who both found that the level of interstellar Call 
absorption is far better correlated with distance for sight-lines 
> lOOpc than for Nal absorption. 

3.3. Column Density Inversion Maps 

Although the plots of equivalent width versus distance can of- 
ten reveal how absorption accumulates over sight-lines (by re- 
vealing gradients in the gas density distribution), the signifi- 
cance of the highest levels of absorption is lost at large dis- 
tances due to line saturation, and similarly information is lost 
on absorption with levels less than that of the lowest measur- 
able limit of equivalent width. However the column density in- 
version method of Vergely et al. (2001) used in Paper 1, which 
derives volume gas densities in 3-D space from the integrated 
sight-line column density values by an iterative inversion pro- 
cess, can reveal masses of gas at any distance provided there 
are sufficient constraining target measurements. 

The inversion algorithm we have presently used to con- 
struct the 3-D gas density maps is derived from a non lin- 
ear least-squares approach to the generaUsed inverse problem 
(Tarantola & ValettelgS^^. It differs from a maximum likeli- 
hood method in that it allows treatment of under-determined 
cases by imposing some additional constraints. For the present 
Nal and Call data we use 'a priori' information on the gas den- 
sity that follows a Gaussian law with non null spatial autocor- 
relation. In other words, the computed gas density distribution 
can be seen as fluctuations around this 'a priori' solution. The 
mathematical specificity is that the solution belongs to the L2 
Hilbert functional space, which allows retrieval of the solution 
everywhere without 'a priori' space discretization. For the case 
of a linear operator linking the unknown parameters (the gas 
densities) and the data (the measured columns) and an assumed 
Gaussian fluctuation around the 'a priori' value, a solution can 
be determined directly in one iteration. However this is not the 
case for our present data (i.e. column densities and distances), 
since in order to force the gas density p to be positive we use a 
change of variable, namely a = \og{p I po) , po being a constant. 
In this case the solution to obtaining the gas density as a func- 
tion of distance is through an iterative method, as described in 
Tarantola & Valette ( 1982) and Vergely et al ( [200TT l. 

The inversion method was first applied to opacities in 3-D 
space by Vergely et al. (1200 11 1. Here we use the same algo- 
rithm, but with the following constraints that are specific to our 
present set of data: (i) the correlation between the gas density 
at point X and the density at point X' is of the form exp(-|X - 
X' p / ^^), with a correlation length, of 20 pc. This distance 
corresponds to the maximal resolution allowed by the present 
data set and is roughly the average distance between two target 
stars. This is an improvement upon our previous inversion of 
Nal data which used ^ = 25 pc, the improvement being allowed 



by the present larger data set. It impUes that details smaller than 
this characteristic scale will be smoothed out and appear at least 
20 pc wide. It also means that the computed volume densities 
are not the local values, but are mean values in volumes whose 
size is of the order of 20^^ pc, (ii) the 'a priori' density decreases 
exponentially with the distance from the galactic plane, with a 
scale height of 170 pc for Na and 450 pc for Call. These two 
values were determined from the data sets themselves by fitting 
an exponential law (for the density) to the measured columns, 
(iii) the 'a priori' error on the opacity, a = \og{p I po), follows a 
Gaussian law with a standard deviation a{a) = 2.5. This allows 
high contrasts for the density, i.e. very dense clouds or close to 
null values, (iv) the errors on the data follow a particular law 
chosen to minimise the effects of potential outliers within the 
data set, as originally devised by Vergely et al. 019981 ), namely 

/(c(*) - co(*)) = TT ■ ait) ■ coshi^-^^^^^-^^) (1) 

where f(c(i)-Co(i)) is the probability for a specific column 
density c(i) towards target i to depart from the central data value 
co(i). Developments show that for close values this function 
tends to a classical Gaussian (a L2 norm), while for values far 
above or below the central one it tends to vary like exp(-(c(i)- 
Co(i))/cr(i)), i.e. a LI norm. 

(v) Relative errors on the columns are taken to be 25% for 
all targets, and those errors are combined with the individual 
errors on the Hipparcos parallaxes. 

Finally, we want to emphasize that within our data set of 
gas densities there are some locations for which there is no 
information that can change the volume density from the initial 
'a priori' value. For these locations the gas density value is kept 
unchanged. We have inluded such information in our maps (see 
below). 

The resultant accuracy of the inversion method has been 
evaluated by Vergely (,1998) through extended tests using syn- 
thetic cloud distributions. The comparison between these ini- 
tial synthetic distributions and the resultant inverted densities 
shows that there are optimal values for both the correlation 
length (also called the smoothing parameter, x) and the statis- 
tical error on the density (i.e. the allowed density fluctuation 
amplitudes). The correlation length is clearly linked to the av- 
erage distance between the target stars and the optimal selected 
value essentially defines the resulting placement accuracy of 
gas clouds using the inversion method. In our case this distance 
is of the order of 35 parsecs for the entire target sample, but is 
less at smaller distances from the Sun. This is because nearby 
targets are more numerous, and thus more information can be 
obtained within the first 150 parsecs. Hence we have selected 
the smaller value of 20 pc for the correlation length in our new 
maps. However, the standard deviation for density fluctuations 
must be non-negligible, since we know in advance that there 
are dense areas and voids between the many groups of dense 
clouds. 

In order for the reader to appreciate the typical accuracy 
of the inversion method for the placement of interstellar gas 
clouds, in Figure 11 we show an example of a fractal distri- 
bution of measurements of interstellar reddening in the galac- 
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tic plane and the corresponding inversion of cloud densities 
computed for 2500 sight-lines randomly distributed in distance 
and directions around the galactic plane. These data are taken 
from Figures 5.10 and 5.11 of Vergely ( |1998l l. with a correla- 
tion length of 15 pc and measurement errors being simulated as 
'noise in the reddening and distance data. Vergely (1998) found 
that the broad structures of reddening (which are equivalent to 
gas clouds) were well located after the inversion, but clouds (or 
voids) smaller than 10-20 pc were included in broader cloud 
structures and thus smeared out. Finally, it must also be men- 
tioned that since the inversion method is statistical in nature, it 
is impossible to prevent some bad restitution for one or more 
cloud structures. This may happen when the combination of the 
actual distribution in one particular galactic area and the target 
sampling in that area are highly non-uniform (i.e. when the re- 
gion or galactic direction is sparsely sampled). 

The end-product of the inversion method, when applied to 
our present measurements (and associated errors) of column 
density and distance for each of the sight-lines sampled, is a 
3-D data cube of Nal (or Call) gas density as a function of dis- 
tance and galactic longitude and latitude. Our new Nal volume 
density data cube has a 60% increase in the number of targets 
over that of Paper I, and the Call maps are the first ones of their 
type ever to be presented. In Paper 1 the correlation length for 
the Nal data was ^ 25pc, whereas the gas distribution shown in 
our new maps of Nal and Call can be considered representative 
at the 20pc sampling scale. Thus, as mentioned above, small- 
scale absorption features (i.e. small cloudlets of low density 
gas) could easily be missed in our final smoothed density maps. 
For a few cases very high column density values of log N(Nal) 
and log N(Call) >> 13.5 cm^^ are reported by other authors 
in Table 1. In some cases it was evident that the absorption 
profiles had not been fitted with sufficient cloud components 
and we noted that such very high density values tended to bias 
and unrealistically distort the 3-D mapping program (with as- 
sociated points of inflection) when low density sight-lines are 
spatially close to these apparently very dense regions. In order 
to smooth out this effect, but still retain the presence of the high 
density regions in our maps, we determined that by generally 
limiting the maximum values of log N(Nal) and log N(Call) to 
13.5 cm^^ this distortion effect was minimized in the resultant 
plots. 

Since we have created 3-D data cubes for both Nal and Call 
gas density out to a distance of 300pc from the Sun, in principal 
we can present maps of the spatial density distribution of either 
ion in any desired galactic projection. However, for brevity, we 
present maps for the 3 main projections (galactic, meridian and 
rotational planes) for both Nal and Call ions. It is envisaged 
that the entire data-base of absorption measurements, together 
with the associated 3-D gas density cubes for Nal and Call will 
be made available to the community on the world-wide web by 
the end of 2009. 

4. Nal absorption within 300pc 

In the following sub-sections we present maps of the 3-D dis- 
tribution of Nal absorption within 300pc of the Sun, formed 
from three main galactic projections through the previously de- 



scribed 3-D Nal gas density cube. Each of the maps represents 
the spatial distribution of the volume density of Nal (n^va/) as a 
function of distance and thus is a 3-D visualization of the data 
previously discussed in Section 3.2. Each map shows the re- 
gions where very low values of neutral gas density are located 
(white shading) together with regions where statistically sig- 
nificant amounts of neutral gas have been detected (grey/black 
shading). Areas with insufficient information to perform an in- 
version of the columns are indicated by a matrix of dots in 
the maps. They correspond to areas within which the density 
is still the 'a priori' density. In cases where some information 
is present but the number of constraints from the target stars 
makes the location of the cloud, or of the cavity, very impre- 
cise, the inversion process may create 'finger-like' extensions 
that are easily recognisable in the maps. The presence of such 
features can be interpreted as indications of the existence of a 
gas cloud along the 'finger' direction, but the location of the 
accumulation of gas is thus uncertain. 

Stars used to create each map are shown as triangles with 
their size being proportional to their distance above (triangle 
vertex pointing up) or below (vertex pointing down) the galac- 
tic plane. Superposed on the maps are contours of iso volume 
density with values of log iiNai = -9.5, -9.1, -8.5 and -7.8 cm^'^ 
shown respectively in yellow, green, turquoise and blue lines. 
As an approximate guide, these average values of neutral gas 
volume density can be converted into hydrogen density values 
at a particular distance using the relationship given in Ferlet et 
al. (1985). Finally, we further remind the reader not to over- 
interpret every absorption feature that appears in these maps 
due to the reasons listed previously. 

4. 1. Nal: A view from above the galactic plane 

The greatest number of sources (i.e. early-type stars) sampled 
by this survey lie close to the galactic plane and hence the gas 
density maps derived for both Nal and Call for this particular 
galactic projection are the most accurate. In Figure 12 we show 
the spatial distribution of Nal density to 300pc which clearly 
reveals a large and highly irregularly shaped volume of space 
surrounding the Sun to a distance of ~ 80pc in most directions 
that is devoid of dense condensations of cold and neutral gas 
clouds. This region, of course, can be identified as the Local 
Cavity (LC). We note that the LC is more elongated in galactic 
quadrants 3 and 4 compared with the other two quadrants, with 
quadrant 1 having the smallest volume of low-density neutral 
gas. Our new map also reveals for the first time that the LC (in 
the galactic plane) can essentially be considered as two rarefied 
'sub-cavities' that abut each other along a thin interstellar gas 
filament of low column density ((log N(Nal) ^ 10.4 cm~^) at 
a distance of ^15pc in the direction of / ^ 345°. One of these 
'sub-cavities' is mostly contained within galactic quadrants 2 
and 3 (and to a far lesser extent within quadrant 1), with the ma- 
jority of the other low density 'sub-cavity' mostly being con- 
tained within quadrant 4. The appearance of smaller 'cell-like' 
cavity structures bounded by thin gas filaments is also hinted 
at in the equivalent maps of Call absorption and these will be 
discussed further in Section 5.1. 
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The relatively dense neutral gas wall that surrounds the LC 
in the galactic plane has a typical depth of 50 -80pc in most di- 
rections, but it is quite evident from the map that the LC void is 
not fully enclosed by a continuous boundary of cold and neu- 
tral gas. Instead, the LC can be best described as being 'porous' 
in the sense that it appears to be linked to several adjacent inter- 
stellar cavities through narrow gaps in the surrounding neutral 
gas wall. These adjacent interstellar cavities include the Loop 
I superbubble, the /3 CMa interstellar tunnel and the Pleiades 
bubble, all of which were revealed in the preliminary maps pre- 
sented in Paper 1 . Our new map now reveals that the tunnel of 
low interstellar gas density towards the direction of the star j3 
CMa ( I = 226°) is bifurcated at a distance of ~160pc, which 
is also the distance to /3 CMa itself. The tunnel's most pro- 
nounced region of low neutral density extends to at least 250pc 
in the direction of I ^ 260°, which provides a pathway to the 
GSH 238h-00h-09 supershell which may have been created by 
the star cluster Collinder 121 (Heiles il99 8). The pathway to 
the Loop I region (I ^ 345°) is now revealed to be far nar- 
rower than that originally shown in Paper 1 , with only one low 
density passageway leading towards the Sco-Cen OB associa- 
tion. The present map shows this narrow entrance to consist of 
several fragmented clouds located over the distance range 90 - 
120pc, in agreement with previous UV observations of this re- 
gion by Welsh & Lallement ( 120051 ). Our new map shows three 
other regions of low neutral gas density (interstellar bubbles?) 
that exist beyond the neutral boundary to the LC at distances > 
200pc. These are located in the directions of I = 190° (Taurus), 
210° (Orion/Eridanus) and 285° (Carina). 

Figure 12 also reveals a previously unknown gap in the neu- 
tral boundary to the LC in the galactic plane that lies between / 
= 70° to 80°. This low density extension of the LC lies in the di- 
rection of Cygnus and may be linked to part of the Loop II SNR 
which lies at a distance of llOpc (Berkhuiisen 119731) . We 
note that in the direction of this newly discovered gap in the LC 
wall there is a large patch of emission revealed in both the 3/4 
keV and 1/4 keV maps of the SXRB (Snowden et al [T997b . The 
presence of such (low neutral gas density) gaps and openings 
in the wall of the LC should have a profound effect on the po- 
tential observability of more distant sources of soft X-ray dif- 
fuse emission. Although this effect is most pronounced when 
viewing the halo through the openings of the Local Chimney 
at high galactic latitudes (see Section 4.2), the gaps may also 
provide a low absorption path for emission from soft X-rays 
generated in nearby hot superbubbles. For example it is well- 
known that there is a dipole effect in the distribution of the soft 
X-ray background emission, with the emission being warmer, 
of greater intensity and with a larger spatial coverage in direc- 
tions towards the galactic center than along sight-lines towards 
the galactic anti-center at mid-plane latitudes (Snowden et al 
12000) . This effect may now be explained through inspection of 
Figure 12 which shows that the local void is linked to the Loop I 
superbubble (which is a known source of million degree soft X- 
ray emission) through a gap in the surrounding boundary wall 
in the direction of the galactic center. 

The closest accumulations of neutral gas in the galactic 
plane (with log N(Nal) ^ 10.8 cm^^) are located in the di- 
rections of (i) I ~ 85° at a distance of ^ 38pc towards the star 



HD 184006, (ii) in the general direction of the galactic center 
at a distance of ^ 50pc (towards HD 159170), and (iii) towards 
several stars with distances in the 60 - 80pc distance range in 
the direction of I ^285°. All 3 of these low neutral density gas 
clouds are extensions of the surrounding denser neutral wall 
that protrude into the central regions of the LC. They appear to 
form the filamentary gas boundaries to the inner cell structures 
whose possible presence within the larger LC was alluded to 
previously. 

In addition, in the direction of / ^^195 ° at a distance of 
^ 95pc we detect another small, more dense gas cloud. The 
presence of two small gaps in the neutral wall to the LC in this 
general direction at / = 165° and / = 190° suggest that the iso- 
lated dense cloud is most probably part of the broken boundary 
wall to the LC. Finally, we note that although our sampling 
of targets with distances in the 200 - 400pc range is far from 
complete (as indicated by regions of small dots in the map), 
there do appear to be galactic directions in the galactic plane 
in which the neutral gas density does not increase with the dis- 
tance sampled. For example, in the direction of / = 315° we 
have an appreciable number of targets sampled with distances 
> 200pc which show no additional Nal absorption beyond that 
distance. This suggests that there may be other large regions 
of low neutral gas density (i.e. superbubbles) that remain to be 
revealed by future more extensive absorption measurements in 
the interstellar region between 200 and 400pc from the Sun. 

4.2. Nal absorption viewed in ttie meridian plane 

In Figure 13 we show the spatial distribution of Nal absorption 
to 300pc as projected in the meridian plane with the galaxy es- 
sentially being viewed in a side-on manner. The plotting sym- 
bols are the same as those for Figure 12. The cold and dense 
Nal absorption associated with gas in the galactic plane can be 
seen as a ~200pc (vertically) thick gaseous bar that stretches 
from left to right across Figure 1 3 . The low neutral density void 
in the central region of this plot is the LC, which in this projec- 
tion can be identified with the open-ended Local Chimney that 
is tilted at an angle of ^ 35° from vertical. This new plot con- 
firms the recent studies of Crawford et al. (2002) and Welsh 
et al. (2004), in which no continuous neutral boundary to the 
LC can be found at high latitudes in either galactic hemisphere. 
This map also confirms the general lack of neutral gas for sight- 
lines that extend to at least '^175pc into the inner halo regions. 
We also note that the rarefied Local Chimney region fragments 
into finger-like extensions that reach into the inner halo in both 
hemispheres. Filaments of neutral HI gas have been widely ob- 
served in the general ISM and their appearance within the local 
ISM suggests that they may be linked to the origin of the LC 
in which an explosive event may have 'cleared out' any neutral 
gas that may have once resided close to the Sun and ejected it 
into the overlying galactic halo. 

Figure 13 shows four major condensations of cold neutral 
gas lying within the confines of the LC in this galactic projec- 
tion. They appear to lie along a fragmented bar of dense gas 
that lies at an angle of ~45° to the galactic plane that stretches 
from 200pc below to 30pc above the plane. The nearest of these 
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clouds lies at a distance of ~ 50pc at a galactic latitude of 6 
30° towards HD 145570 and it is an extension of the more dis- 
tant dense wall to the LC in that direction. The three remaining 
neutral clouds appear as a string of gaseous 'blobs' in the direc- 
tion of (I ~ 165°, b ~ -75°) at distances of 60pc (towards HD 
9672), 1 lOpc and 180pc. The cloud of cold neutral gas at 1 lOpc 
can be associated with the translucent molecular cloud G192- 
67 (Grant & Burrows I19991 I. This string of 3 nearby clouds 
would be an ideal candidate for future soft X-ray shadow obser- 
vations that potentially could resolve the outstanding issue of 
the degree of pervasiveness of the purported million degree gas 
that may be present within the Local Cavity (Welsh and Shelton 
[2009). Finally, Figure 13 also shows a thin and elongated cloud 
of neutral gas at a distance of ~ 200pc in the galactic direction 
of (/ ~ 335°, b ^ -30°), whose existence is revealed for the 
first time. 

4.3. Nal absorption viewed in ttie rotational plane 

Figure 14 is perhaps the most changed of all three neutral 
density maps when compared to the previous presentations of 
Paper 1, due to the small number of target stars available along 
this plane in the previous work. The low neutral density rar- 
efied LC region extends to ^ 60pc in most directions, and 
the surrounding fragmented wall of much denser neutral gas 
is generally not encountered until a distance of lOOpc. As 
also revealed in the meridian projection of Figure 13, there is 
no continuous neutral boundary to the LC at high latitudes in 
both galactic hemispheres. In the rotational projection we see 
that the Local Chimney has only one major narrow opening to- 
wards the north galactic pole (tilted at b ~H-75°), whereas the 
opening into the galactic halo in the southern hemisphere is 
of far larger dimensions. Figure 14 also shows narrow exten- 
sions of the rarefied inner LC to distances as far as ~ 200pc 
in several directions. Their finger-like shape means their dis- 
tances are not well constrained, but it is interesting to try and 
relate them to known distant features. The most prominent of 
these cavity extensions are seen towards the southern galactic 
pole and also towards / ^ 250° close to the galactic plane. The 
latter feature is the pathway to the GSH 238h-00h-09 supershell 
discussed previously in Section 4. 1 . 

Figure 14 suggests that there are no large accumulations of 
neutral gas contained within the confines of the rarefied LC in 
this projection. The nearest edge of a Nal cloud is at a distance 
of ^ 40pc in the direction of I ^ 90° near to the galactic plane 
(towards the stars HD 184006 and HD 186882). This cloud ap- 
pears to be an extension of the denser neutral wall to the LC 
in this galactic direction. We further note the presence of two 
denser Nal clouds at distances of ~ 85pc and 1 lOpc at galactic 
latitude 6 ^75°. Both clouds are also probably fragments of 
the surrounding disk gas that may have been disrupted by the 
expansion of the original Local Bubble into the overlying less 
dense galactic halo. 

One particularly interesting feature shown in this map is 
the circular region of low neutral gas density at a distance of 
~140pc that lies 25° above the galactic plane in the general 
direction of (/ ^ 90°). The low level of absorption along this 



sight-line was first reported by Lilienthal et al. (I1991I I and we 
now confirm the presence of a large interstellar cavity with sim- 
ilar dimensions to that of the Loop I superbubble feature. 

For completeness we note that Meyer et al. ( 120061 ) have re- 
ported the detection of a very cold (T ^ 20K) neutral Nal cloud 
located at a distance of 45 - 55pc in the direction (I ~ 230°, 
b H-45°). Unfortunately the galactic position of this cloud is 
such that it does not appear in any of the 3 galactic projections 
we have currently presented. However, this is an important de- 
tection since the cloud's presence within the Local Cavity has 
been posed as a challenge to current theories concerning the 
survival of cold gas clouds that may be immersed in a purport- 
edly hot (million K) surrounding plasma. 

5. Call absorption within 300pc 

In the following sub-sections we now present three maps of the 
spatial distribution of Call absorption within 300pc of the Sun, 
formed from a 3-D Call volume density cube derived in a sim- 
ilar manner to that described for the Nal data. The new maps, 
derived with a correlation length of 20pc, are shown in Figures 
15-17 and have the same plotting symbols as those presented 
for Nal, except that the Call volume density iso-contours (yel- 
low, green, turquoise and blue) correspond to values of log 
ncaii = -9.9, -9.5, -8.9 and -8.2 cm^^ . 

5. 1. Call: A view from above the galactic plane 

As was the case for the Nal data, this map has the largest num- 
ber of sight-line measurements of Call and thus is the most 
accurate of the 3 different galactic projections all presented 
here for the first time. The overall appearance of Figure 15 
shares a gross similarity with the equivalent plot for Nal ab- 
sorption in that it reveals a central region of mainly low den- 
sity Call gas that is surrounded by a highly fragmented wall 
of denser partially ionized gas clouds. However, on closer ex- 
amination of Figure 15 we see that (in common with the Nal 
data) the LC is not one large volume of rarefied low neutral 
density gas, but in fact it can be better described as a collection 
of several smaller low density cavities that are each surrounded 
by thin gas filaments of slightly higher density neutral and/or 
partially ionized material. Although some of these filamentary 
structures may be artifacts caused by uncertainties in the in- 
version mapping method when only a few target measurements 
are present along an given sight-line, this is not the case for all 
of the diffuse gas filaments that appear to be present in Figure 
15. Clearly, although more measurements will be required to 
confirm the reality of a collection of cell-like cavity structures 
that may be present within the LC, such a scenario would have 
some intriguing astrophysical ramifications. For example, how 
could such small interstellar cell-like cavities form within the 
Local Cavity? One possibility is that each of the rarefied sub- 
cavities could be a small ionized stellar HII region (with a neu- 
tral/partially ionized filamentary gas envelope) formed by the 
stellar wind action of nearby early B-type stars (and to a lesser 
extent by local hot white dwarfs). We note that there are at least 
1 1 stars of spectral type earlier than B3V located within lOOpc 
of the Sun whose combined stellar UV flux may be responsible 
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for photo-ionization of the surrounding low density medium. 
This scenario would provide a natural explanation for the ob- 
served ionization state of the local interstellar medium as mea- 
sured over many local sight-lines by Lehner et al. (2003). We 
note that previous estimates of the ionization sources of the lo- 
cal ISM have assumed the presence of a hot million degree gas 
in addition to photoionization from the early-type star e CMa 
in order to explain the physical state of the local interstellar 
cloud (Slavin & Frisch 2008). Unfortunately a detailed calcu- 
lation and comparison of the local ionization field is beyond the 
scope of this present work, and we urge theorists in this field 
to test whether the picture of a Local Cavity with no million 
degree gas (Welsh & Shelton I2009I I. but instead containing a 
collection of small ionized HII region cavities can reproduce 
the UV absorption line strengths and velocity structure widely 
observed along local sight-lines by Lehner et al. (1 20031 ) and 
Welsh & Lallement ([2005]) . 

In contrast with the equivalent absorption map of Nal, the 
central region of the LC has at least six condensations of ap- 
preciable Call absorption (i.e. clouds of of partially ionized 
gas) lying within 80pc of the Sun. The nearest region of Call 
absorption is the complex of local interstellar clouds that ex- 
tend to ^ 15pc mostly in directions contained within galac- 
tic quadrants 1 and 4. Redfield & Linsky (I2008I I. using both 
ground-based and UV absorption observations recorded to- 
wards nearby stars, have presented more detailed maps of each 
of the several clouds that contribute to this local cloud complex 
and argue that the Sun is located in a transition zone between 
two of these clouds. Our results are in quantitative agreement 
with their results, with the majority of this local Call absorp- 
tion being attributed to the presence of the local interstellar 
cloud, the G-cloud and the Blue cloud (Lallement et al. 1995 
Redfield & Linsky 2008). The second nearest region of appre- 
ciable Call absorption is located at a distance of ^50pc in the 
general direction of (I ^80°) towards the stars HD 186882 and 
HD 193369. This feature is also seen in the corresponding Nal 
absorption map as a very low neutral density finger-like gas 
cloud that looks like an extension of a larger and denser inter- 
stellar cloud located at ^ lOOpc in the same galactic direction. 
The third nearest region of Call absorption is located at ^ 55pc 
in the direction of (/ -325°) towards the stai's HD 134481, HD 
134482 and HD 141194. It has no apparent Nal counterpart. 
The fourth nearest region of Call absorption is located at a dis- 
tance of — 60pc in the direction of{l^ 20°) and this may either 
be an extension of the local interstellar cloud complex in this 
direction or a protrusion from the surrounding wall to the LC. 
Finally, there are two regions of appreciable Call absorption 
located at about 70pc towards / ~ 265° and I 110°. While 
the former cloud has a slightly more distant Nal counterpart, 
none is clearly seen for the second regions. This gas may be 
strongly photo-ionized, as is the case for more distant gas seen 
in this direction (see next paragraph). 

Figure 15 shows that the dense Call (and Nal) cloud struc- 
tures that surround the LC are accompanied by many narrow 
pathways (tunnels) of lower Call gas density that extend into 
surrounding interstellar regions. Many of the dense regions of 
Nal absorption that surround the LC close to the galactic plane 
(as shown in Figure 11) are also seen as high density Call- 



bearing clouds at the same physical locations. In addition, most 
of the interstellar tunnels revealed by the paucity of Nal ab- 
sorption in Figure 12 (and discussed in Section 4.1) are also 
observed as low density regions in the Call absorption map of 
Figure 15. The major exception to this is the low neutral den- 
sity CMa tunnel (Z = 226°) that is so prominent in the Nal maps. 
Figure 15 reveals this direction to possess a high level of Call 
gas density, which is to be expected if this tunnel of interstel- 
lar gas is photo-ionized by the two early-type stars of e and /3 
CMa (Vallerga 1998). We also note that the narrow interstellar 
tunnel of low neutral gas density leading into the more distant 
cavity at (/ — 345°) is significantly more obscured in the Call 
maps. In contrast, the region lying beyond 150pc in this direc- 
tion is devoid of both neutral and partially ionized gas, as one 
might expect if the cavity (which we associated with the Loop 
I superbubble) is hot and highly ionized. 

The main conclusion from Figure 15 is that there are a sig- 
nificant number of partially ionized diffuse gas clouds that are 
present within the LC. Their ubiquity, current state of ioniza- 
tion and gas temperature have yet to be fully accounted by any 
model of the local ISM (Slavin & Frisch EUQHI VallergaT998l 
Bruhweiler & Cheng 1988) . and thus we urge a re-assement of 
some of these models in the Ught of the recent findings of the 
lack of a local million degree gas and the possible existence of 
cell-hke cavity structures within the LC. 

5.2. Call absorption viewed in the meridian plane 

In Figure 16 we show the spatial distribution of Call absorp- 
tion within 300pc as projected in the galactic meridian plane. 
The central region of low Call absorption density essentially 
mimics the Local Chimney feature shown in Figure 13 for Nal 
absorption, but the corresponding Call Chimney has a narrower 
diameter of — lOOpc and reveals several small cell-like struc- 
tures within the rarefied central region. The LC (when viewed 
in Call) is also seen to be open-ended towards both high and 
low galactic latitudes, with a larger (fragmented) opening to- 
wards the southern hemisphere, in concert with the correspond- 
ing map of Nal neutral gas absorption. The central low den- 
sity volume of the LC is surrounded in the galactic plane by 
dense partially ionized Call gas clouds mostly located within 
the galactic disk. Many of these clouds are also seen in the 
equivalent Nal absorption map of Figure 13. However, we note 
that the narrowing of the Chimney in Figure 16 is mainly due 
to the presence of dense Call clouds located in galactic quad- 
rants 2 and 3 whose full extent are not traced by Nal absorp- 
tion. For example, the string of 3 gaseous Nal 'blobs' seen in 
the direction of (? - 165°, 6 - -75°) in Figure 12 have a Call 
counterpart that has the appearance of a large hook-like cloud 
that comes as close as — 40pc to the Sun. 

There are 3 major condensations of partially ionized Call 
gas lying within lOOpc of the Sun in this projection. The near- 
est of these is the local cloud complex (d < 15pc) that partially 
surrounds the Sun and has been discussed in the previous sec- 
tion. We note that, as suggested by the galactic plane projection 
map, this cloud appears to be physically linked to the wall of 
denser gas that surrounds the LC region in the direction of the 
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galactic center. The second closest cloud of high Call gas den- 
sity lies at a distance of ^70pc in the direction of the north 
galactic pole. This cloud has no apparent Nal counterpart. The 
third closest dense cloud of Call is located at a distance of 
85pc with a galactic latitude, b ~ -50°, and again it has no ap- 
parent counterpart in the equivalent Nal absorption map. 

In the direction of the south galactic pole we find a region of 
high Call absorption at a distance of ~ 130pc, whose position 
is coincident with that of a finger-like extension of low den- 
sity neutral gas seen in Figure 13. On the other hand, the large 
circular Call cloud in the general direction of I ^ 180° and dis- 
tance ~ 150pc with a galactic latitude of 6 +45° has no Nal 
counterpart. Due to the lack of suitable targets with distances > 
200pc at high galactic latitudes, we are unable to trace the dis- 
tribution of partially ionized gas further into the lower galactic 
halo. Clearly, additional observations of fainter and more dis- 
tant early-type stars will be required in order to accomplish this 
task. 

5.3. Call absorption viewed in the rotational plane 

Figure 17 shows the spatial distribution of Call absorption 
within 300pc in the galactic rotational plane projection. We 
observe a large circular volume of very low Call absorption 
density that surrounds the Sun to a distance of ^ SOpc in most 
directions. This low density cavity has an SOpc wide opening 
into the overlying galactic halo at high positive latitudes, but 
there is only a narrow tunnel of low Call gas density that leads 
into the inner halo towards the southern galactic pole. Within 
the central low density region there are only two clouds of high 
density Call absorption, both of which seem to be physically 
part of a bar of partially ionized gas that includes the local in- 
terstellar cloud (within which the Sun is partially immersed). 
One of the high density regions of partially ionized gas is lo- 
cated at ~ SOpc in the galactic plane in the direction of I ^ 90°, 
and can be identified with the spatially coincident gas cloud 
detected in Nal and discussed in Section 4.3. The other local 
high density Call cloud is in the general direction of / ^ 270° 
at a distance of ^ 35pc, located just below the Sun in Figure 
17. This cloud complex can be identified with absorption due 
to both the interstellar G-cloud and Auriga cloud as listed by 
Redfield & Linsky (120081) . We also note that the bar of Call gas 
within the LC seems to be linked to the surrounding boundary 
wall (at least in the direction of / ^ 90°). Its physical shape 
and dimensions are suggestive of a cell-like structure, in agree- 
ment with the local web of low-density interstellar sub-cavities 
discussed in Section 5.1. 

The fragmented wall of high Call gas density that mostly 
surrounds the central LC region is generally matched with 
spatially coincident regions of high Nal neutral gas density. 
However, this one-to-one matching is not always observed. 
For example, the tunnel of low neutral density gas that opens 
into the galactic halo towards the north galactic pole shown in 
Figure 14 is seen to be spatially coincident with a high density 
region of Call absorption. Similarly, the large high density Call 
cloud at 1 lOpc seen towards the south galactic pole is matched 
by a similarly large region of very low Nal gas density. This 



pattern of Nal and Call absorption behavior is reflected in the 
derived Nal/Call column density ratio, whose major variations 
within 300pc are discussed in the following section. 

6. Variation of thie Nal/Call ratio 

The column density ratio of N(Nal)/N(Call) is a well-known 
diagnostic of the physical conditions in the diffuse interstellar 
gas, since Ca is more sensitive than Na to the balance between 
adsorption on, and desorption from, interstellar grains (Barlow 
119781) . In cold (T ^ 30K) and dense gas clouds, in which most 
of the gas-phase Ca is depleted onto grains, the Nal/Call ratio is 
> 100. However, in the warmer (T lOOOK) and lower density 
ISM where much of the Ca remains in the gas phase, ratios of 
< 1.0 are commonly found (Hobbs'1975! Centurion & Vladilo 
119911 Bertin at al. 119931 ). For sight-line distances < 30pc the 
Nal/Call column density ratio for the warm (T ^ 7000K) local 
interstellar gas clouds is ~ 0.2 (Bertin et al. il993) . Equilibrium 
equations show that at 7000K the Nal/Call ratio is independent 
of electron density and as a consequence this ratio becomes a 
tracer of the level of calcium elemental depletion for the local 
gas. 

Beyond the LC and in the general ISM the Nal/Call ratio is 
thought to be a function of gas cloud velocity, such that at ve- 
locities 30 km s^^ interstellar dust grains may be destroyed 
by shocks and the Ca is liberated into the gas phase (Routly & 
Spitzer [T952l Siluk & Silk [1974^. Thus, as has been argued by 
Crawford et al. ( 120021 ). a low value of the Nal/Call ratio could 
therefore be due to either the presence of warm and partially 
ionized gas and/or the presence of interstellar shocks. Since we 
do not expect, nor have observed, interstellar gas with veloci- 
ties ^ 30 km s~^ in the local ISM then for distances ^ lOOpc 
any variation in the observed Nal/Call ratio is most probably 
due to temperature/ionization effects. From our sample of stars 
that possess measured values of both Nal and Call column den- 
sities we see that there is a range in Nal/Call ratio values of 
0.03 to 100, with the vast majority of sight-lines having ratio 
values < 5.0. Since the column density values listed in Table 
1 are total integrated values measured over the entire length of 
each sight-line, the derived Nal/Call ratio value is a sight-line 
average value. Thus, variations in the ratio for low column den- 
sity cloud components that may be present along a particular 
sight-line will be masked by the column density contribution 
from any high column density (cold and dense) regions. One 
example of this is shown in Figure 2 for HD 301 12, whose to- 
tal Nal/Call column density ratio is 2.1 whereas the V 5 km 
s^^ cloud component has a ratio of only 1.15. Thus, we defer 
a detailed study of the spatial variation of the Nal/Call ratio 
in the local ISM to the forthcoming paper by Lallement et al. 
(12010 ) that will focus on the cloud component velocity struc- 
ture of these sight-lines. 

In Figure 1 8 we show plots of the column density ratio of 
Nal/Call as a function of distance to 400pc for sight-lines with 
latitudes near the galactic plane {b ±20°) and split by longitude 
into 4 galactic quadrants. We see that for distances < 80pc the 
majority of measured values of the Nal/Call ratio lie in the 0. 1 
to 1.0 range, in general agreement with the low value for this 
ratio found by Bertin et al. ( 119931 ) for stars within the Local 
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Cavity. For more distant sight-lines the Nal/Call is more vari- 
able, but in most directions the ratio value falls in the range 
0.5 to 20. This is in agreement with with previous work on the 
Nal/Call ratio in the general ISM by Welty et al. (1 19961 1. 

The highest values of the Nal/Call ratio near the galactic 
plane are found in quadrant 2 in the direction of 1 ^ 150° to- 
wards the Taurus dark clouds (i.e. along sight-lines towards 
HD 23180 (55.0), HD 23552 (72.4), HD 24398 (85.1) and HD 
25642 (102.3). There appears to be qualitative agreement in the 
distribution of the Nal/Call with distance in all 4 galactic quad- 
rants, with a smaller range of ratio values (typically 0. 1 to 5) 
present in quadrant 3. This is most probably due to the pres- 
ence of significant amounts of partially ionized Call gas in the 
200pc long low neutral density interstellar tunnel seen towards 
the stars (3 and e CMa which shows an ionization gradient in 
this galactic direction (Wolff et al. |1999l l. 

7. Conclusion 

We have presented new high spectral resolution absorption 
observations of the interstellar Nal D-lines (5890A) recorded 
along 482 sight-lines, together with 807 new measurements of 
the Call K-line (3933A). When combined with existing data 
in the literature, we have compiled a catalog of absorption 
measurements towards a total of 1857 early-type stars located 
within 800pc of the Sun (1678 stars measured at Nal and 1267 
stars measured at Call). By measuring both the Nal and Call 
interstellar line equivalent widths and then subsequently fitting 
the normalized absorption profiles with a line velocity, doppler 
width and column density, we have been able to characterize 
the neutral and partially ionized properties of the interstellar 
medium along each sight-line. A plot of the equivalent width 
of Nal D2-line versus distance reveals a central region of very 
low neutral gas absorption (Wa(D2) < 5mA) out to a distance 
of 80pc in most galactic directions. This is the well-known 
Local Cavity region, which has a wall of neutral Nal gas (of 
Wa(D2) > 50mA) beyond 80pc in most directions. In contrast, 
a similar plot for the equivalent width of Call shows no sharply 
increasing absorption at 80pc, but instead we observe a slowly 
increasing value of Call equivalent width with increasing sight- 
line distance sampled. 

Low values for the volume density of Nal (n„a/ <10^^ 
cm~^) are generally found within 50pc of the Sun (i.e. within 
the Local Cavity), whereas values in the range 10~^ > njva/ > 
10^^° cm^"^ are found for sight-lines with distance > 300pc. 
Both high and low values of the volume density of Call (ncaii) 
are found for sight-lines < 30pc, dependent on whether local 
gas cloudlets are encountered. For distances > lOOpc a value 
of ncaii ^ 10^^ cm~^ is typical for most sight-lines, suggest- 
ing that the density distribution of Call gas absorption is fairly 
uniform throughout the general ISM. 

Each of the normalized absorption profiles of the Nal and 
Call lines were then subsequently fit with models of line com- 
ponent velocity, doppler width and column density, with the to- 
tal integrated values of column density for each ion along each 
sight-line being presented in Table 1 . These new data were then 
used in conjunction with previous absorption measurements 
pubhshed in the literature to construct 3-D maps of the spa- 



tial distribution of neutral Nal and partially ionizes Call gas 
density. This was achieved from the inversion of each column- 
density determination using a method based on the work of 
Vergely et al. (2001 ). The Nal maps extend and improve upon 
the accuracy of similar maps constructed by Lallement et al. 
(12003), whereas the Call maps are the first of their kind to be 
published. 

The general appearance of the maps of the 3-D spatial dis- 
tribution of Nal gas density are similar to those of Call absorp- 
tion, in that we find a ~ 80pc diameter central region of low 
absorption (i.e. the Local Cavity) surrounded by a wall of high 
absorption density. In contrast to the preliminary maps of Nal 
absorption presented in Paper 1, the surrounding wall of Nal 
(and Call) absorption is now found to be highly fragmented 
with several new interstellar tunnels of low gas density ema- 
nating from the central LC region. This fragmentation could 
be linked to the origin of the LC, which is generally thought 
to have been formed by a SN explosion over several million 
years ago. For most sight-lines many regions of high Nal ab- 
sorption are spatially matched by regions of high Call absorp- 
tion. However, there are several important differences within 
the details of these maps. Firstly, the central region of the Local 
Cavity contains several clouds of high Call gas density, in con- 
trast to the general absence of neutral Nal gas clouds found 
within the rarefied cavity. The observed pattern of local Call 
absorption appears to be in the form of a collection of sev- 
eral very low density interstellar cells that are surrounded by 
partially ionized Call filaments. This collection of low den- 
sity cavities might be caused by the combined action of stellar 
winds from nearby early B-type stars that could be responsible 
for the photoionization of the local interstellar gas. Secondly, 
although there are low gas density pathways that lead from the 
LC into the overlying galactic halo, the low density neutral Nal 
gas openings are of far larger dimensions than their low density 
Call equivalents. Thirdly, there are several regions of high Nal 
absorption that do not have high Call absorption counterparts, 
and vice-versa. This is reflected in the value of the integrated 
column density ratio of Nal/Call for these clouds, which is pre- 
sumably influenced by the ambient level of ionization in these 
regions. 

Plots of the integrated column density ratio of Nal/Call as 
a function of distance for sight-lines near to the galactic plane 
have values that lie in the range 0. 1 to 1 .0 for sight-lines with 
distances < 80pc. However, ratio values of between 0.5 and 
20 are typical for more distant sight-lines. The highest values 
of the Nal/Call ratio are found towards I ^ 150° in the direc- 
tion of the Taurus dark clouds, and ratio values in the narrower 
range of 0.1 to 5 are found in galactic quadrant 3 due to the 
presence of the low neutral density /3 CMa interstellar tunnel. 
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Fig. 1. Examples of interstellar Call-K line absorption pro- 
files together with best-fit models (thick line) recorded with 
the Aurelia spectrograph on the 1.52m telescope at the 
Observatoire de Haute Provence. Dotted lines are the model 
components prior to instrumental convolution. The best-fit 
component values of velocity (V), doppler-width (b) and col- 
umn density (N) are listed in the box beneath each profile. 
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Fig. 2. Examples of interstellar Nal D2 and Dl and Call-K 
line absorption profiles together with best-fit models (thick 
line) recorded with the Hamilton spectrograph on the 0.9 m 
CAT telescope at the Lick Observatory of the University of 
California. Dotted lines are the model components prior to in- 
strumental convolution. The best-fit component values of ve- 
locity (V), doppler-width (b) and column density (N) are listed 
in the box beneath each profile. 
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Fig. 3. Examples of interstellar Nal D2 and Dl and Call-K line 
absorption profiles together with best-fit models (thick line) 
recorded with the Hercules spectrograph on the 1.0 m telescope 
at the Mt. John Observatory in New Zealand. Dotted lines are 
the model components prior to instrumental convolution. The 
best-fit component values of velocity (V), doppler-width (b) 
and column density (N) are listed in the box beneath each pro- 
file. The short vertical lines above some profiles indicate the 
position of the absorption component's velocity. 



Welsh & Lallement: Nal and Call maps 



1.0S 




-AO -20 20 40 
Heliocentric vriocity (km/s) 



1.05 
S 1.00 
S 0.9S 
3 0.90 
0.85 ^ 
0.80 



t ' ■ ■ 1 ' < < 1 ■ ' * 1 ■ ' 






Dl 


- 

HD 480 





-40 -20 20 40 
Heliocentric velocity (km/s) 



V,= 
v,= 



+1.9 b, =3.4 logN: = 10.7li0.15 
+SJ b. =3.4 logN.= I0.49±0.1i 




-40 -20 20 40 
Heliocwitrtc veiocity (km/s) 



V,--K>.S b, -2.9 logNi - 1 1.08±0.05 
Vj = -i-7.0 bn = 2.9 logN2= 1 1.15+fl.06 



1.0 
' 0,8 



I 0.4 
I 0.2 
0.0 





\ 


V : 




HD 1 34796 



-40 -ZO 20 40 
Heliocentric velocity (km/s) 




•40 -20 
Heliocentric velocity (km/s) 




-40 -20 20 40 
Melwcentrtc velocity (km/s) 



V,- 


• 17.7 


b, 


-3.6 


logN, 


- II.23±0.IS 


V:- 


•8,2 




-3.S 


k>gN: 


- 10.93±0.12 


V, - 


■1,4 


b, 


- 1.5 


logN, 


• 1 1.84*0.1,1 


v., 


IH.I 


h 


2.7 


lug N. 


! 1 I7i« 15 



Vi 


-22.0 


b, 2,9 


IukN, - I0.9JtO.I: 


V. • 


-1,2 


b; - 2.9 


lOSN;- ll.49±O IO 


V, 




b, 1.S 


litgNi - 12 tlKt.2n 




-40 -20 20 40 
Heliocentnc velocity (km/s> 



1.04 - 

I t.oo 

c 

1 0.96 

■g 

= 0.92 



' ' 1 ' ■ ' 1 ' ' ■ 1 


■ ' ' 1 ' ' ' 1 ' ' 










■ Dl 

.1 .... 1 


HD 48261 ■ 



-40 -20 20 40 
Heliocentric velocity (km/s) 



V, 


= +12.6 


b, 


= 3.2 


log Ni 


= 10.59+0 12 


V, 


-+2i5 




-3.3 


IoeN, 


- 10.68±0.15 


V, 


=+2S.S 


bi 


'13 


log Na 


= 10.1K±0.20 



1.04 



1.00 



0.96 



0.92 



Cal-K 




-40 -20 



HO 48261 



20 40 



Hefcocentric vefocity (km/s) 



V,=+IO,0 


b. 


= 2.9 


log N, 


= 1 0.65+0. 12 


V. =+I7.6 


b; 


= 2.3 


log N. 


- 10.60*0.14 


V3=+25.l 


bj 


= 2.9 


log N, 


= 11.03+0. IS 



Fig. 4. Examples of interstellar Nal D2 and Dl and Call K and 
H- line absorption profiles together with best-fit models (thick 
line) recorded with the FEROS spectrograph on the 2.2 m tele- 
scope at the European Southern Observatory in Chile. Dotted 
lines are the model components prior to instrumental convolu- 
tion. The best-fit component values of velocity (V), doppler- 
width (b) and column density (N) are listed in the box beneath 
each profile. 
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Fig. 5. Examples of interstellar Nal D2 and Dl and Call-K line 
absorption profiles together with best-fit models (thick line) 
recorded with the GIRAFFE spectrograph on the 1.9 m tele- 
scope at the South African Astronomical Observatory. Dotted 
lines are the model components prior to instrumental convolu- 
tion. The best-fit component values of velocity (V), doppler- 
width (b) and column density (N) are listed in the box beneath 
each profile. 
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Fig. 6. Galactic distribution of the interstellar sight-lines sam- 
pled by (a) Nal absorption, and (b) Call absorption. Both plots 
are centered at galactic co-ordinate (0,0). 
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Fig. 7. Plot of the equivalent width (mA) of the interstellar Nal 
D2-line for stars with distances < 400pc. Filled triangles are 
for sight-lines with galactic latitude b > +45°, open squares 
for sight- lines with 6 = to 45°, open circles for sight-lines 
with 6 = to -45° and filled circles for sight-lines with b < 
-45°. Crosses are upper limit values. Note the sharp increase 
in the level of Nal absorption at ~ 80pc, which is due to the 
neutral wall to the Local Cavity. 
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Fig. 8. Plot of the equivalent width (mA) of the interstellar Call 
K-line for stars with distances < 400pc. Symbols are the same 
as those in Figure 7. 
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Fig. 9. Plot of the volume density of interstellar Nal absorp- 
tion for stars with distances < 800pc. Symbols are the same as 
Figure 7. 
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Fig. 10. Plot of the volume density of interstellar Call absorp- 
tion for stars with distances < 800pc. Symbols are the same as 
those in Figure 7. 
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Fig. 11. Comparison of an initial map of the fractal distribu- 
tion of interstellar reddening in the galactic plane (left) and the 
equivalent map of gas column densities (right) derived from 
the inversion method of Vergely (1998). Most large features 
are reproduced, whereas small features with sizes < 15pc are 
smeared by the inversion method. 
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Fig. 12. Plot of 3-D spatial distribution of interstellar Nal ab- 
sorption within 300pc of the Sun as viewed in the galactic plane 
projection. Triangles represent the sight-line positions of stars 
used to produce the map, with the size of the triangle being 
proportional to the derived Nal column density. Stars plotted 
with vertex upwards are located above the galactic plane, ver- 
tex down are below the plane. White to dark shading repre- 
sents low to high values of the Nal volume density(njva/)- The 
corresponding iso-contours (yellow, green, turquoise and blue) 
for log njva/ = -9.5, -9.1, -8.5 and -7.8 cm^'^ are also shown. 
Regions with a matrix of dots represent areas of uncertain neu- 
tral gas density measurement. 



26 Welsh & Lallement: Nal and Call maps 




Fig. 13. Plot of 3-D spatial distribution of interstellar Nal ab- 
sorption within 300pc of the Sun as viewed in the meridian 
plane projection. Plotting symbols are the same as those given 
for Figure 1 1 . 
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Fig. 14. Plot of 3-D spatial distribution of interstellar Nal ab- 
sorption within 300pc of the Sun as viewed in the rotational 
plane projection. Plotting symbols are the same as those given 
for Figure 1 1 . 
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Fig. 15. Plot of 3-D spatial distribution of interstellar Call ab- 
sorption within 300pc of the Sun as viewed in the galactic plane 
projection. Plotting symbols are the same as Figure 11, except 
that the volume density iso-contours (yellow, green, turquoise 
and blue) correspond to log ncaii = -9.9, -9.5, -8.9 and -8.2 
cm~^ 
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Fig. 16. Plot of 3-D spatial distribution of interstellar Call ab- 
sorption within 300pc of the Sun as viewed in the meridian 
plane projection. Plotting symbols are the same as Figure 14. 




Fig. 17. Plot of 3-D spatial distribution of interstellar Call ab- 
sorption within 300pc of the Sun as viewed in the rotational 
plane projection. Plotting symbols are the same as Figure 14. 
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Fig. 18. Plots of the column density ratio of N(NaI)/N(CaII) 
as a function of sight-line distance for galactic quadrants 1-4. 
Black circles are measured values, open circles are upper limits 
and crosses are lower limits. 
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Table 1. Nal and Call Absorption Measurements 



Star//D 


1 


b 


distance 


Wa (D2) 


Wa (D1) 


log N(Nal)roT 


Reference 


Wa (K) 


log N(CaIl)Tf)r 


Reference 








(pcj 


ni(A) 


m(A) 


cm " 


(Nal) 


m(Aj 


cm " 


(Call) 


179029 


0.1 


-20 


318 


124.3 


105.4 


12.37 


ESO-2007 


53 


12.01 


ESO-2007 


157056 


0.5 


6.6 


173 


87.0 




11.98 


(1) 


38.0 


11.67 


(2) 


172016 


0.8 


-12.5 


225 


136.6 


76.20 


11.97 


ESO-2007 


25.7 


11.33 


ESO-2007 


195599 


0.9 


-36.1 


186 


<2.5 




<10.15 


ESO-2007 


3 


10.65 


ESO-2007 


203006 


1.2 


-45.0 


57 










<5 


<10.4 


(3) 


171034 


1.3 


-11.0 


685 


240.0 


175.0 


12.39 


MJO-2007 


83.4 


12.06 


MJO-2007 


152909 


1.6 


14.4 


209 


235.0 


198.5 


12.83 


ESO-2006 


104.9 


12.24 


ESO-2006 


154204 


1.9 


12.4 


122 


179 


142 


12.66 


(4) 








164019 


1.9 


-2.6 


556 










498 


12.99 


(5) 


161756 


2.0 


0.5 


348 


284.8 


236.6 


12.87 


ESO-2006 


99.9 


12.30 


ESO-2006 


145570 


2.5 


28.8 


51 


6.7 


6.1 


10.67 


ESO-2008 


21.9 


11.52 


ESO-2008 


158643 


2.5 


5.3 


131 






11.05 


(6) 


20.3 


11.26 


(7) 


165365 


2.8 


-3.7 


418 


186.4 


121.8 


12.15 


ESO-2008 


58.3 


11.91 


ESO-2008 


180885 


2.8 


-20.7 


565 


109.8 


86.7 


12.16 


MJO-2006 








197630 


2.9 


-38.3 


98 


12.9 


3.3 


10.71 


MJO20-06 


10.3 


11.06 


MJO-2007 


151884 


3.2 


17.4 


268 


372.8 


348.9 


12.89 


ESO-2008 


131.8 


12.38 


ESO-2008 


45607 


3.9 


29.7 


77 


30.0 


15.0 


11.28 


(8) 


17.6 


11.36 


LICK-2007 


141378 


4.1 


37.3 


49 


<3 




<10.3 


(9) 








141569 


4.2 


36.9 


99 






12.66 


(10) 




11.65 


(10) 


188113 


4.2 


-27.8 


446 


199.5 


152.6 


12.50 


ESO-2008 


66.9 


12.0 


ESO-2008 



